****************************************
* Table B6: Democrat Vote Share U.S. House: Closest Distance
******************************************

preserve
use "$input/prec-closest-hur", clear
xtset precid



qui gen treatment = treatment50
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6a.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)   replace dec(3)   
drop cmean coef percent treatment
local foo


qui gen treatment = treatment100
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6a.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)   append dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment200
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6a.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)   append dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment300
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6a.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)   append dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment400
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6a.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)   append dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment500
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6a.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)    dec(3)   append tex
drop cmean coef percent treatment
local foo

use "$input/prec-closest-major-hur", clear

qui gen treatment = treatment50
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6b.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)  replace  dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment100
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6b.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)  append  dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment200
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6b.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)  append  dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment300
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6b.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)  append  dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment400
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6b.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)  append  dec(3)   
drop cmean coef percent treatment
local foo

qui gen treatment = treatment500
qui xtreg demohouseshare_total treatment i.year i.stateid, fe cl(precid)
qui gen coef=_b[treatment]
qui sum demohouseshare_total if treatment==0 & e(sample)
qui gen cmean = r(mean)
qui  gen percent = coef/cmean * 100
local mean1 = round(cmean, 0.001)
local percent1 = round(percent, 0.001) 
outreg2 using "$output/tabb6b.doc", addtext("Control Mean =", "`mean1'", "Effect  =", "`percent1'",)   dec(3)   append tex
drop cmean coef percent treatment
local foo


restore


